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Abstract 

New superintegrable systems are presented which, like the Hydrogen atom, possess a 
dynamical symmetry w.r.t. algebra o(4). One of them simulates a neutral fermion with 
non-trivial dipole moment, interacting with the external e.m. field. This system is pre- 
sented in both non-relativistic and relativistic formulations. Another recently discovered 
system (see arXiv:1208.2886vl) is non-relativistic and includes the minimal and spin-orbit 
interaction with the external electric field. It is shown that all the considered systems 
are shape invariant. Applying this quality, these systems are integrated using the tools of 
SUSY quantum mechanics. 
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1 Introduction 



The Hydrogen atom (HA) is one of the most important systems of quantum mechanics. This is a 
perfect physical model with a nice symmetry. Namely, in addition to its transparent invariance 
with respect to the rotation group, the HA possesses a hidden (Fock [1]) symmetry w.r.t. group 
0(4) whose generators are the orbital momentum and Runge-Lenz vector [2j. Moreover, this 
system is also supersymmetric and can be solved algebraically using tools of SUSY quantum 
mechanics [3]. 

At the best of my knowledge, until the first version of this paper appeared in ArXiv- 
1205.3094vl there were only two known generalizations of the Runge-Lenz vector for 3d QM 
systems with spin. The first of them was discovered by Johnson and Lippman as far back as 
in 1950 [1] as a hidden symmetry of the relativistic Dirac equation with Coulomb potential. A 
contemporary treatment of the Johnson-Lippman constant of motion in arbitrary dimensional 
space is presented in [5]. The other generalization was proposed in 1985 by D'Hoker and L. 
Vinet who proved that the Schrodinger-Pauli equation for electron interacting with the dyon 
field admits a vector integral of motion depending on spin [6]. However, analogy of these 
symmetries with the Runge-Lenz vector is rather poor. Indeed, the Johnson-Lippman integral 
of motion is a scalar which is decoupled to three vector components only in the non-relativistic 
limit. The constant of motion found in |6f includes the same higher derivative terms as the 
Runge-Lenz vector, but it does not generate the dynamical symmetry w.r.t. group 0(4). 

In very recent paper [9] superintegrable QM systems with spin-orbit interaction and second- 
order integrals of motion are classified. One of these systems which includes the generalized 
Runge-Lenz vector is shortly discussed below in Section 6. 

Notice that the spin dependent Runge-Lenz vector was discovered also for for spinning 
Taub-NUT space [7], but it was done on the pseudoclassical level. 

In the present paper a new QM system with spin 1/2, which admits a hidden symmetry 
w.r.t. group 0(4), is discussed. Its integrals of motion are the total orbital momentum and 
the generalized Runge-Lenz vector depending on spin. Like the HA, this system appears to be 
supersymmetric also, which makes it possible to find its exact solutions using the regular SUSY 
approach. 

Mathematically, the models presented in the following are interesting new examples of su- 
perintegrable and supersymmetric systems with spin, which are exactly solvable. These systems 
include shape invariant matrix potentials, which are particular cases of such potentials classified 
recently in [TO] , [H]. In addition, these examples corroborate the conjecture that all maximally 
superintegrable systems are exactly solvable p2] (see [13] , [H] , [15] for discussion) and present 
a new field for studying the relations between supersymmetry and superintegrability. 

Physically, the models discussed in sections 3-5 simulate a neutral fermion with a non- 
trivial dipole moment (e.g., the neutron), interacting with the external field. Moreover, a 
relativistic version of this model is presented also. It is shown that the neutron can be trapped 
by the specific external fields. Potentially, these systems and their exact solutions can have 
a wide spectrum of applications, starting with using them as tutorial examples and ending 
with relevance to the problems of security of nuclear reactors. In any case, the superintegrable 
systems like the HA and isotropic harmonic oscillator are extremely important in physical 
applications. Thus we can hope that, to some degree, it would be the case for the systems 
discussed in the present paper. 
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2 Symmetries of the Hydrogen atom 



Let us remind the main symmetry properties of the HA which will be used as a standard for 
construction of the model with spin. 

The Hamiltonian of the HA looks as follows: 

H=^- q - (1) 

2m x 



where p 2 = p\ + p\ + p§, p\ = — « gf-, % = + x \ + x i' ? > 0. 

Operator ([1]) is manifestly invariant w.r.t. rotation group 0(3) and so it commutes with 
the angular momentum vector: 

L = x x p. (2) 

In addition, it commutes with the Runge-Lenz vector 

R = — (p x L — L x p) + xV (3) 

2m 

where V — — -. 

X 

Components of vectors fl2]) and ([3]) satisfy the following commutation relations 

[L a , Lb] = ie a b c L c , [R a , L^\ = ie a b c R c , 

2i (4) 

[Ra,Rb] — £abcL c H. 

m 

On the set of eigenvectors of Hamiltonian H algebra (jl]) is isomorphic to the Lie algebra of 
group 0(4). 

Thus there are the six integrals of motion for Hamiltonian ([T]) (we do not take into ac- 
count the Hamiltonian itself). However, we can select maximally four such integrals which are 
algebraically independent and include a pair of commuting representatives. This means that 
the HA is a maximally superintegrable system: the number of its algebraically independent 
constants of motion including the Hamiltonian is equal to 2n — 1, i.e., is maximal for a system 
with n = 3 degrees of freedom. 

The other (discrete) symmetry of Hamiltonian (Tj[|) is the space reflection P. Indeed, this 
Hamiltonian is transparently invariant w.r.t. the change x — > — x. 

In addition, operator ([Q) admits a hidden supersymmetry [3j. In other words, its radial 
component is shape invariant with respect to the special Darboux transformation which will 
be specified in Section 4. 



3 Runge-Lenz vector for fermions 

Operators (J2]) and ([3j) by construction correspond to a spinless system. A straightforward way 
to generalize them for the case of a fermion system is to change orbital momentum L by the 
total angular momentum: 

L ->• J = L + S. 

1 ~ (5) 

R R = (p x J - J x p) + xV 
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where S = |cr is the spin vector and cr is the matrix vector whose components are Pauli 
matrices. Potential V should commute with J and is in general spin dependent. 
Operators (j^J) should satisfy relations (jlj) where L a — > J a , R a — )■ R a : 

[J a, Jb] '^abcJci [Ra, Jb\ ^abcRc, 

~ . 2i (6) 

[Ra, Rb] — ZabcJcH. 

m 

and commute with a Hamiltonian which we will search in the form 

H = L + *- < 7 > 

It can be verified by direct calculation that up to a constant multiplier a we have the unique 
choice for V: 

V = a°-±. (8) 

x z 

Thus it is Hamiltonian with "matrix Coulomb potential" OH]) which admits the integrals 
of motion This Hamiltonian is not invariant w.r.t. space inversion x — > — x since cr is a 
pseudovector. However the space inversion is an admissible transformation provided the wave 
function ip cotransforms in a non-standard way: 

V>(*,x)-r£V>(*,-x) (9) 

where B = * s ^ e Biedenharn operator [T6]. More exactly, the following condition is 

satisfied: 

OT(p,x)5 = #(-p,-x). 

A natural question whether this Hamiltonian succeeds another symmetry of (jTJ), i.e., the 
shape invariance, is discussed in the following section. 

4 Supersymmetry and exact solutions 

Consider the eigenvalue problem for Hamiltonian ([7]): 

£. + a ^-) * = s*. do) 

Introducing rescalled independent variables r = 2max it is possible to rewrite (flUI) in a more 
compact form: 

-A + ^)^ = ^ (11) 



where e = • 

Taking into account the invariance of equation (flOj) w.r.t. the rotation group, it is convenient 
to rewrite it in the spherical coordinates and expand ip y i a spherical spinors Qjj-\ )K ((p, 6): 

^ = ~ S $i^( T )®ij-*JjP* d )- ( 12 ) 



Here j = ~, |, k = —j, + ... j, and A = ±| are quantum numbers labeling 

eigenvalues of the commuting operators J 2 , L 2 and J 3 : 

L 2 Q j7j -x, K = (j - A)(j? - A + l)n 3 -j_A,K, 
The explicit form of the spherical spinors is given by the following formula [17] : 




a., i = 2 2 , a,, i 
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where ^±1^1 are spherical functions. 

Substituting ( fl2|) into ( TTTj) we come to the following equations 

U 1 



n ^ =\-9^+Vi) ^ = ^ ( 13 ) 



where 



and V ? - is the matrix potential: 



Vi=(j(j + l)+\-<T 3 (j+lj) 



--en-. (15) 



Equations (IT51) appear to be supersymmetric since Vj belongs to the list of shape invariant 

potentials classified in [TO], see Eq. (5.11) for k = |, /i = j there. Indeed, Hamiltonian can 
be factorized as 

ft,- = a+Oj + c, (16) 

where 

a ' = Tr + w '- «t—h + w * c ^-TUTTf (17) 

and Wj is the matrix superpotential 

"H^-H^^V 1 - (18) 

Moreover, Hamiltonians Tij and "Hj+i satisfy the following intertwining relations 

H ja + = atHj+i. (19) 
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Thus Eq. (JTHjl can be easily solved using tools of SUSY quantum mechanics. In fact it has 
been already done in [TO]. The ground state $° fc = column (-0°, £°) should solve the first order 
equation a^- k = 0. Thus, using definitions (117)1 we obtain: 



*!.* = * .., 3 „ 2 °„ +1) < I (20) 



2(i+i) 



where and Kq are the modified Bessel functions, Ck are arbitrary constants. The corre- 
sponding eigenvalue Eq is equal to Cj, i.e., 

1 



4(j + l) 2 ' 



The solution $™ fc corresponding to n th exited state and the corresponding eigenvalue e n are: 



®lk ~ a t a t+i ' ' ' °i+n-i^i+n,fe! £ n - 4 / J - + n + 1 \ 2 - ( 21 ) 
The related energy value in ( ITUj) is given by the following equation: 

where 

JV = n + j + l, n = 0,1,2,... (23) 

The energy levels (I22p depend on the main quantum number iV ( I2"3"j) which can take the 
same value for different pairs of n and j. Namely, for a fixed A" there are A" — | of such pairs. 
In addition, the quantum number k which labels the eigenvectors in (1131) is not present in (T22|) 
and ( 123]) . Thus like in the HA problem the energy levels are highly degenerated, and this 
degeneration is caused by the hidden symmetry w.r.t. group 0(4). 

To end this section we remind that the radial component of the HA Hamiltonian (P) is 
shape invariant also and can be factorized like f lT6|) where j — y I, W — > 2 {i+i) ^ anc ^ ^ * s 
the quantum number labeling eigenvalues of angular momentum L. In other words, our model 
succeeds both hidden symmetry and supersymmetry of the HA. 



5 Relativistic system 

Let us show that eigenvalue problem (fTTj) admits a relativistic formulation. To demonstrate 
that we start with the Dirac equation for a neutral particle having non-trivial dipole moments: 

( 7 % - m - aa^F^ + aj^Fj V = 0. (24) 



In addition to the standard Pauli term aa^ v F^ v Eq. (1241) includes additional term aj^j^F^ 
with an external pseudovector field F M . For convenience the following realization of Dirac 
matrices will be used: 

70 = U o) ' 75 = U ~o) ' 7a = 0? -Ij (25) 
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where / and are the 2x2 unit and zero matrices correspondingly, a = 1, 2, 3. 
Choosing in 



F ab = F = 0, F a = ^F 0a = -f (26) 
a x 2 

and representing ^ as: 

* = ex P (-^x ) (^\ (27) 

with two- component functions if) and £, we come to the following system: 

(itr ■ p - m)ip + E£ = 0, (28) 

E - 2a ^r) $ ~ (™> + icr ■ P)C = 0. (29) 

Solving (|28j) for ^ and substituting this solution into (|29|) we come to the eigenvalue problem 
(TTTT) where r = 2aEx and e = ■ Thus, in accordance with ( 12T|) the relativistic energies 

are discrete also and can be given by the following formula: 

E2 = Ss (30) 

with jV being a natural number which can be represented in the form (|23j) . The corresponding 
(non-normalized) wave functions are given by Eq. (j27]) where ^(x) is defined in (fl2|) . (fl4l) . 
pi. pi]) and ^(x) can be found from ([28]). 

Notice that for small coupling constants a the positive energy values (13T)j) are reduced to 
the nonrelativistic form ( 12 2 p up to the constant internal energy term m and terms of order a 4 . 



6 Superintegrable and supersymmetric system with spin- 
orbit interaction 

In very recent preprint [12] one more QM model with a hidden symmetry w.r.t. group 0(4) is 
presented. The Hamiltonian of this model looks as follows: 



v 2 - l 

H = -f- + V - A, 31 

2m 8x 2 



where 



^=w( CT - L + 1 >-i- (32 > 

In contrast with ([7]) Hamiltonian fl3Tl) does not include the dipole interaction term but 
involves a scalar potential together with the spin-orbit interaction. Nevertheless, it also com- 
mutes with the total orbital momentum and generalized Runge-Lenz vector §5§ where V is 
given by equation f[3"2"j) (in paper [H] another representation for R is used which is equivalent 
to ([5]), (|32j) ). This result was predictable since Hamiltonian (13T|) can be reduced to the direct 
sum of two Hamiltonians of HA via the gauge transformation [8] . 
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It can be proven by direct verification, that operators (jSJ), ( 132]) and ( 13T|) satisfy commutation 
relations (EJ) and so they generate a hidden symmetry with respect to group 0(4). Thus the 
QM system with Hamiltonian (131 p admits a generalized Fock symmetry. 

Hamiltonian f l3~T|) succeeds also one more symmetry of the systems considered in the above, 
i.e., the shape invariance. Indeed, starting with f l3~T]) and repeating all steps following equation 
ffTTj) . one comes to radial equation (TT3]) where 



V, = ^- \. (33) 

Here and in the next equation each term in the r.h.s. is multiplied by a 2 x 2 identity matrix. 

Like (fl5]l . potential f )33|) is shape invariant. The related superpotential is given by the 
following formula: 

1 7 + 1 

W = - r + J —— (34) 

20 + 1) r K J 

while eigenvalues of Hamiltonian (13 ip are given by equation (j22p where n and j are positive 
integers and half integers correspondingly. The corresponding wave function is easily calculated 
in complete analogy with section 4. It is given by equation ( 112"]) where ij}j,x,k V^A,fc anc ^ 

^A, fe = ca^' +1 exp (-|) L? +1 (»). (35) 

Here y = n ^f +1 , L% +1 (y) are Laguerre polynomials and c\k are constants satisfying the nor- 
malizing condition Yl\k c ^k c *\k = 1- 

Effective potential ( 13"3"|) and superpotential (|34p are similar to ones for the HA. However, 
there are some differences, namely: 

• Potential fl33l includes quantum number j which takes half-integer values while the po- 
tential of the HA depends on orbital quantum number I which is integer; 



• In contrast with the HA, potential fl33l is a matrix. More exactly, it is a direct sum of 
two scalar potentials. That leads to the additional two-fold degeneration of the energy 
spectrum of Hamiltonian f[3"T]) . This degeneration is caused by the additional integral of 
motion 

C = - J • R = ^ (36) 

whose eigenvalues are not included into the spectrum formula (!22j) . 

• Matrix (136]) is proportional to the Casimir operator C = J • R of group 0(4), and its 
eigenvalues are equal to ±|. For the HA this Casimir operator is trivial. 

• The additional symmetry operator (I3"6l) extends the number of algebraically independent 
constants of motion (including the Hamiltonian) to 6. This number is maximal for a 
3d system with the additional (spin) degree of freedom, thus the discussed system is 
maximally superintegrable. 

The matrix fl3"6"|) and many other integrals of motion for Hamiltonian fl3"T|) were represented 
in paper [Bj. All of them are algebraic functions of the basic symmetry operators fl5]) where V 
is the matrix given in (|32i) . 
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7 Discussion 



The model Hamiltonian (JJJ represents a new exactly solvable QM system which admits extended 
symmetries. Like the HA this system admits the hidden (Fock) symmetry w.r.t. group 0(4) 
whose generators are the total angular momentum and Runge-Lenz vector. In addition (and 
again like the HA) this system is supersymmetric and can be easily solved using tools of SUSY 
QM, see Section 4. Moreover, this is a fermionic system with spin ~ while the non-relativistic 
model of HA ignores the spin of electron. 

A new feature of the supersymmetry of Hamiltonian ([7]) in comparison with the HA is that it 
is realized in terms of matrix superpotentials. Such (one dimensional) superpotentials had been 
classified in papers [10J and and it was very inspiring for us to search for multidimensional 
matrix systems which can be solved using their matrix supersymmetry in separated variables. 
As a result Hamiltonian ([7]) has been discovered. The other results of such search are presented 
in papers [H] and [i~9] . 

Operator ([7J) can be interpreted as the Hamiltonian of a neutral particle with spin 1/2 
which has a non-trivial dipole moment. Matrix potential = otcr ■ E represents a Pauli 

type interaction with the external vector field 

e~4- ( 37 ) 

x z 

The latter can be interpreted as the electric field. Moreover, such field can be realized exper- 
imentally at least on finite interval a < x < b, a > 0, see, e.g., problem 1018 in [20]. Notice 
that functions ( 137]) also solve nonlinear Maxwell equations including additional vector field 
[2"T] . They also solve equations of axion electrodynamics [22]. The same is true for fields (T26]) 
involved into relativistic equation (T2"4"|) . 

Thus the QM systems whose Hamiltonians are given by equations ([7]) and ( 13Tj) posses all 
symmetries and supersymmetries admitted by the HA. However these symmetries are general- 
ized by introducing the spin. In addition, there exist a relativistic counterpart of system (TlO|) 
given by Eqs. ( 124")) and ( 126]) which is exactly solvable too. 

The presented non relativistic system f llOp can be treated as a 3d generalization of the 
solvable planar system proposed by Pronko and Stroganov [23] ■ It would be interesting to 
study its possible generalizations in multidimensional spaces. One more interesting task is to 
construct solvable systems with Fock symmetry for particles with arbitrary spin. This work is 
in progress. 

It is well known that there exists a tight coupling between the hidden symmetry of the HA 
and its supersymmetry [21], [25]. The same is true for the Pronko-Stroganov problem [26J. 
A contemporary discussion of relations between supersymmetry and superintegrability can be 
found in [27] . 

All systems discussed in the present paper are both supersymmetric and maximally su- 
perintegrable. Surely, these combinations of symmetries are not accidental, and it would be 
interesting to extend the results [23], [26J to the 3d matrix systems presented here. 
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